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Abstract. Critical behaviour of a fluid (binary mixture or liquid crystal), subjected to strongly 
anisotropic turbulent mixing, is studied by means of the field theoretic renormalization group. As 
a simplified model, relaxational stochastic dynamics of a non-conserved scalar order parameter, 
coupled to a random velocity field with prescribed statistics, is considered. The velocity is taken 
Gaussian, white in time, with correlation function of the form oc S(t — i ' ) / 1 k jj d +£ , where 
is the component of the wave vector, perpendicular to the distinguished direction ("direction of 
the flow") — the d-dimcnsional generalization of the ensemble introduced by Avellaneda and 
Majda [Commun. Math. Phys. 131 381] within the context of passive scalar advection. It is 
shown that, depending on the relation between the exponent £ and the space dimensionality d, 
the system exhibits various types of large-scale self-similar behaviour, associated with different 
infrared attractive fixed points of the renormalization-group equations. In addition to well known 
asymptotic regimes (model A of equilibrium critical dynamics and passively advected scalar 
with no self- interaction), existence of a new, non-equilibrium and strongly anisotropic, type of 
critical behaviour (universality class) is established, and the corresponding critical dimensions 
are calculated to the second order of the double expansion in £ and e = 4 — d (two-loop 
approximation). The most realistic values of the model parameters (for example, d = 3 and the 
Kolmogorov exponent £ = 4/3) belong to this class. The scaling behaviour appears anisotropic 
in the sense that the critical dimensions related to the directions parallel and perpendicular 
to the flow are essentially different. The results are in qualitative agreement with the results, 
obtained in experiments and simulations of fluid systems subjected to various kinds of regular 
and chaotic anisotropic flows. 
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1. Introduction 

Various systems of very different physical nature (ferromagnets and antiferromagnets, gas-vapour 
systems, binary liquid mixtures and alloys) reveal interesting singular behaviour when undergoing 
continuous (second-order) phase transition (that is, in the vicinity of their critical points). 
Specific heat, susceptibility, spontaneous magnetization etc exhibit singular self-similar (power-law) 
behaviour, whose quantitative characteristics (critical dimensions and scaling functions) depend 
only on few global characteristics of the system (like symmetry or space dimensionality). This 
universality is related to the existence in such systems of a wide range of strongly coupled degrees 
of freedom: it produces a kind of collective behaviour in which numerous irrelevant details of a 
specific system are wiped away. This classical subject is exposed in the monographs pQ and the 
literature cited therein. 

Consistent qualitative and quantitative description of the equilibrium critical behaviour was 
achieved within the framework of the renormalization group (RG) . In the RG approach, possible 
types of critical regimes (universality classes) are associated with infrared (IR) attractive fixed 
points of renormalizable field theoretic models. Most typical phase transitions belong to the 
universality class of the 0(./V)-symmetric (p 4 model of an iV-component scalar order parameter 
(Landau-Ginzburg Hamiltonian). Universal characteristics of the critical behaviour depend only 
on N and the space dimensionality d and can be calculated in the form of the expansion in e = 4 — d 
or within other systematic perturbation schemes; see the monographs 0E] and the literature cited 
therein. 

Dynamical critical behaviour (critical singularities of relaxation and correlation times, various 
kinetic and transport coefficients etc) appears richer, less universal and is comparatively less 
understood. Different nature of the order parameter (conserved or non-conserved), inclusion of 
"secondary" slow modes (densities of entropy or energy) and interaction with hydrodynamical 
degrees of freedom produce different types of critical dynamics for the same static model EI El ■ 
The reliable sets of second-order (two-loop) results were only recently fixed, and many important 
questions remain open; see the recent review paper 5 for discussion and bibliography. 

It has long been realized that the behaviour of a real system near its critical point is 
extremely sensitive to external disturbances, geometry of the experimental setup, gravity, presence 
of impurities and so on. "Ideal" equilibrium critical behaviour of an infinite system can be obscured 
by limited accuracy of measuring the temperature, finite-size effects, finite time of evolution 
(ageing) and so on. What is more, some disturbances (randomly distributed impurities in magnets 
and turbulent mixing of fluid systems) can produce completely new types of critical behaviour 
with rich and rather exotic properties (e.g., expansion in y/s rather than in e); see [HllZIIHl- Over 
the past three decades, considerable attention has been attracted by the effects of various kinds of 
imposed flows (laminar shear flows, turbulent stirring and other types of deterministic or chaotic 
flows) on the behaviour of critical fluids, e.g., of binary liquid mixtures near the consolution point; 
see the papers |7]-[101 an d references therein. This problem is closely related to another interesting 
issue: the effects of imposed flows on the dynamics of phase ordering — the growth of order through 
domain coarsening (spinodal decomposition), when a system is quenched from its high-temperature 
homogeneous phase into the low-temperature multi-phase coexistence region; see e.g. JE]-E3 and 
the literature cited therein. 

Although very different, many of such systems exhibit a common interesting feature: existence 
of new non-equilibrium stationary states with (arguably) self-similar statistical properties and new 
sets of scaling exponents. Emergence of nonequilibrium steady states appears rather a generic 
and robust phenomenon, being observed in experiments, simulations and analytical treatments of 
critical liquids with passive and active order parameters subjected to laminar or turbulent flows, 
various kinds of regular and chaotic synthetic velocity ensembles, cellular or shear flows, and so 
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on [7]-|2U|. In the presence of a distinguished direction, scaling behaviour of such systems appears 
strongly anisotropic, with different critical dimensions corresponding to different spatial directions 

G5HIEI 

The aim of the theory is to establish existence of such regimes on the basis of microscopic 
dynamic models, to classify corresponding universality classes, to calculate their scaling dimensions 
within consistent approximations or regular perturbation schemes, to investigate their universality, 
dependence on the model parameters, and so on. In this paper, we will focus on the anisotropic 
turbulent mixing of critical fluids, because most real flows are strongly anisotropic, this anisotropy 
persists in the asymptotic critical regime and leads to new interesting effects. 

The full-scale model of a critical fluid subjected to a strongly anisotropic turbulent stirring 
must deal with a conserved (binary mixtures) or non-conserved (liquid crystals) order parameter 
with mutual coupling with the velocity field, governed by nonlinear dynamic equations (e.g. 
stochastic Navier-Stokes equation with an external random stirring force), and the anisotropy 
is introduced by the initial and/or boundary conditions. However, even for the equilibrium and 
isotropic case (model H in the traditional classification introduced in 0]), the consistent RG analysis 
of such problem appears a most difficult task, and has only recently been completed (see discussion 
and references in |S] and sections 5.23-5.25 in book ; 3,), while theoretical description of fully 
developed turbulence on the basis of dynamic equations remains, in many respects, essentially an 
open problem. 

In the present paper, we apply the field theoretic RG to a simplified "minimal" model of a 
stirred critical fluid, which nevertheless appears rather nontrivial and captures the main property 
of the problem: existence of a new, non-equilibrium and strongly anisotropic, universality class of 
scaling behaviour. Namely, we consider a purely relaxational dynamics of a non-conserved passive 
scalar order parameter (model A in terminology of 0]) coupled to the random velocity field with 
prescribed Gaussian statistics. 

Recently, the models involving passive (no feedback on the velocity) linear (no self-interaction) 
scalar fields advected by such "synthetic" velocity ensembles attracted enormous attention among 
the "turbulent community" because of the insight they offer into the origin of intermittency and 
anomalous scaling in the real fluid turbulence; see the review paper |21j and references therein. 
In spite of their relative simplicity, such models reproduce many of the anomalous features of 
genuine turbulent heat or mass transport observed in experiments. Most popular is the Kazantsev- 
Kraichnan ensemble with the velocity correlation function of the form (vv) oc 6(t — t')k~ d ~^. 
Vanishing of the correlation time is necessary to ensure Galilean symmetry of the problem, while a 
power-law dependence on the wave number k mimics real self-similar properties of fully developed 
turbulence. For a conserved order parameter, it can be shown that the nonlinearity in the Navier- 
Stokes equation and a finite correlation time are indeed IR irrelevant (in the sense of Wilson) in 
the analysis of critical behaviour; see also the discussion in [Sj- 

In the RG approach to the Kraichnan model and its descendants, reviewed in |22| . the exponent 
£ plays the role of a formal RG expansion parameter, analogous in this respect to the conventional 
e = 4 - d. 

Synthetic ensembles also allow one to easily introduce anisotropy, compressibility etc, and to 
study their effects on the behaviour of the scalar field. In this paper we employ the d-dimensional 
generalization of a strongly anisotropic ensemble introduced in |25 in connection with the passive 
linear problem: the velocity field is oriented along a chosen direction n and its correlation function 
depends only on the coordinates perpendicular to n; see also Refs. |2^1l23] . 

In experiments, critical fluids with a non-conserved order parameter can be realized in twisted 
nematic liquid crystals; see the discussion and references in |17H18| . In a wider context, such model 
can be viewed as a model system for studying generic nonequilibrium phase transitions. Recently, 
significant progress has been achieved in classifying IR scaling behaviour of such phenomena, 
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including driven diffusive systems, diffusion-limited reactions, growth, ageing and percolation 
processes, and so on; see e.g. Refs. and references therein. Being analytically tractable, 

our model can serve as a possible testing ground in studying such scaling regimes and their 
universality within controlled approximations or a regular perturbative scheme. 

Earlier, the field theoretic RG was applied to the stirred critical fluid in a number of studies. 
New types of critical behaviour were identified for model B in a Gaussian velocity ensemble 
IB] and generalized model A with inclusion of a large-scale stirring force and the velocity field 
governed by the stochastic Navier-Stokes equation ^] EOj , but only purely isotropic situations 
were considered. The RG ideas were also applied to the problem of phase separation and domain 
growth below the critical temperature; see e.g. 123 GUI and references therein. By contrast with 
critical phenomena, the RG approach to such problems suffers from the lack of an (obvious) small 
parameter (like e = 4 — d or £ in our case) and should involve numerical (Monte Carlo) simulations 
|22 or additional phenomenological hypotheses . 

The paper is organized as follows. In section |3 we give detailed description of the model, 
present its field theoretic formulation and the corresponding diagrammatic technique. In section|31 
we analyze canonical dimensions and ultraviolet (UV) divergences of the model. We show 
that, after an appropriate extension, the model becomes multiplicatively renormalizable and 
present the corresponding renormalized action functional. We also show that, for the extended 
model, independent canonical dimensions should be introduced for the directions parallel and 
perpendicular to the flow. In section 0] we derive the differential RG equations, introduce the 
RG functions (j3 functions and anomalous dimensions 7) and give the corresponding one-loop and 
two-loop expressions for the case of an iV-component order parameter. In section we analyze 
possible scaling regimes of the model, associated with the fixed points of the RG equations, and 
identify their ranges of stability in the e— £ plane. Three fixed points correspond to known regimes: 
free (Gaussian) regime, linear passive scalar and equilibrium model A. The fourth fixed point 
corresponds to a new, non-equilibrium and strongly anisotropic, universality class. For the most 
realistic values of the model parameters (namely, d = 3 and < £ < 2; see section 0) it is the 
latter point that is IR attractive and governs the large-scale, long-time behaviour of the system. 
The corresponding critical dimensions depend on e and £ and can be systematically calculated as 
double series in those parameters; the explicit second-order results are presented in sectional The 
scaling regime appears strongly anisotropic in the sense that the critical dimensions related to the 
directions parallel and perpendicular to the flow are different. Section[7|is reserved for discussion, 
comparison to the existing experimental and theoretical results and the conclusions. 

In appendix A we explore consequences of the Galilean symmetry for the renormalization of 
our model. The main points concerning the calculation of the renormalization constants and RG 
functions are briefly discussed in appendix B. 

2. The model. Field theoretic formulation 

Relaxational dynamics of a non-conserved scalar order parameter fix) with x = {t, x} is described 
by a stochastic differential equation 

aodMx) = -6H(<p)/5<p(x) + f(x), (2.1) 

where d t = d/dt, <7o = 1/T is the reciprocal of the (constant) kinetic coefficient T > and f(x) 
is a Gaussian random noise with zero mean and the pair correlation function 

D f (x - x>) = (f(x)f(x')) = 2a S(t - f ) ^(x - x'), (2.2) 
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d being the dimensionality of the x space. Near the critical point, the Hamiltonian Ti.(^p) is taken 
in the Landau-Ginzburg form 



%)= / dx|-^(x)aV(x) + ^^ 2 (x) + ^^(x)| 



(2.3) 



where d{ = d/dx{ is the spatial derivative, d 2 = didi is the Laplacian, tq oc (T — T c ) measures 
deviation from the critical temperature and Ao > is the coupling constant; after the functional 
differentiation in l|2.1|l one has to replace </?(x) — > (p(x). The model H2.lfl ~ H2.3fl * s referred to as 
model A 0]; its critical behaviour is very well understood 

Coupling with the velocity field Vi{x) is introduced by the replacement 

d t -> V* = fit + Vid h (2.4) 

where V* is the Lagrangian (Galilean covariant) derivative. 

Let n be a unit constant vector that determines distinguished direction ("direction of the 
flow"). Then any vector can be decomposed into the components perpendicular and parallel to 
the flow, for example, x = x^ + nxu with xj_ • n = 0. The velocity field will be taken in the form 

v = u + m>(i,x ± ), (2.5) 

where u is a constant vector parallel to n and v(t, xjj is a scalar function independent of X|| . Then 
the incompressibility condition is automatically satisfied: 

d i v i = d l] v(t,x 1 _)=0. (2.6) 

From now on, we set u = (the general case u ^ leads to no serious alterations in the RG 
analysis and will be briefly discussed in the end of section[7J). For v(t, xj_) we assume a Gaussian 
distribution with zero mean and the pair correlation function of the form: 

(«(t, xjXf, x'J) - S(t - t') / exp {ik ■ (x - x')} D v (k) = 



(2TT) d 

= 5(t-t')J-^ ZI exp{ik ± -(x x -x> ± )}D v (k ± ), fcj. = |kj.| (2.7) 

with the scalar coefficient functions of the form 

D v (k) = 2ir8(k\ l ) D v (k±), D v (k ± ) = D kl d+1 ^. (2.8) 

Here Dq > is a constant amplitude factor and £ an arbitrary exponent, which (along with 
the conventional e = 4 — d) will play the part of a formal RG expansion parameter. The IR 
regularization in l|2.7fl is provided by the cutoff k± > m (by dimension, tq oc m 2 ). [Precise form of 
the IR regularization is inessential; sharp cutoff is the most convenient choice from the calculational 
viewpoints. Another possibility is to replace k± — > \/k\ + m 2 in If2.8fl .] The natural interval for 
the exponent is < £ < 2, when the so-called "effective eddy diffusivity" 

V(r±) = / (2^T i 1 - ex P ( ik ^ ■ r ^)i (2- 9 ) 

has a finite limit for m — > 0; it includes the most realistic Kolmogorov value £ — 4/3. The exponent 
£ can also be viewed as a kind of Holder exponent, which measures "roughness" of the velocity 
field [23; the "Batchelor limit" £ — > 2 corresponds to smooth velocity. 

In order to ensure multiplicative renormalizability of the model, it is necessary to split the 
Laplacian in 1)2.311 into the parallel and perpendicular parts d 2 — > d\ + uqO 2 by introducing a new 
parameter uq > (in the anisotropic case, these two terms will be renormalized in a different way). 
Thus equation l|2.1fi becomes 

o-qVMx) = dl<p(x) + u d 2 tp(x) - T ifi(x) - X ifi 3 (x)/6 + f(x); (2.10) 
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this completes formulation of the model. 

Interpretation of the splitting of the Laplacian term in (|2.10|) can be twofold. On the one hand, 
the fluctuation models of the type 1)2-1)) and ()2.3)1 are phenomenological and, by construction, they 
must require all the IR relevant terms allowed by symmetry. The fact that the splitting is required 
by the renormalization procedure means that it is not forbidden by dimensionality or symmetry 
considerations and, therefore, it is natural to include the general value uq ^ 1 to the model from 
the very beginning. On the other hand, one can insist on studying the original model with uo = 1 
and SO(d) covariant Laplacian term, although that symmetry is broken to SO(d — 1) by the 
interaction with the anisotropic velocity ensemble. Then the extension of the model to the case 
uo ^ 1 can be viewed as a purely technical trick which is only needed to ensure the multiplicative 
renormalizability and to derive the RG equations. The latter should be then solved with the 
special initial data corresponding to Uo = 1 (in renormalized variables this anyway will correspond 
to general initial data with a / 1). Since the IR attractive fixed point of the RG equations is 
unique for any given choice of the parameters e and £ (see section [SJ) , the resulting IR behaviour 
will be the same as for the case of the extended model with general uq ^ 1. 

According to the general theorem [3^ (see also the monographs HIED, our stochastic problem 
is equivalent to the field theoretic model of the extended set of fields <& = {if', ip, v} with action 
functional 

S($) = (T {ip') 2 + if' [-a Vtip + diip + uodftp-T (p- Xoip 3 /6] + S»(v). (2.11) 

The first few terms represent the De Dominicis-Janssen action functional for the stochastic problem 
H2.ll) . 1)2. 2|) at fixed v; it involves auxiliary scalar response field ip'(x). All the required integrations 
over x — {t, x} and summations over the vector indices are implied, for example, 

<p'dj_(p= dt dxtp' (x)d^(p(x). 



It is worth noting that, owing to transversality of the velocity field 1)2.6)1 . the derivative in the 
coupling term in l|2.11|l can also be moved onto the field ip' using integration by parts: 

ip' '(vidi)tp — dt dxip' (x)v(t,x±)dnip(x) = 



dt J dx(d\i<p'(x))v(t,x.±)<p(x). (2.12) 

The last term in (|2.11l) corresponds to the Gaussian averaging over v with correlator (|2.7I) 
and has the form 

5„(v) = ~ J dt J dx ± dxlw(t,x ± )5- 1 (x ± - x' ± )v(t,x' ± ), (2.13) 

where 

5- 1 (r ± )oc£> - 1 rf- d) - S (2.14) 

is the kernel of the inverse linear operation D" 1 for the correlation function D v in 1)2. 8JI . 

This formulation means that statistical averages of random quantities in the original stochastic 
problem coincide with the Green functions of the field theoretic model with action H2.ll)) . given by 
functional averages with the weight exp £(<!?) (see equation HA. 2)1 in the appendix A). This allows 
one to apply the field theoretic renormalization theory and renormalization group to our stochastic 
problem. The model 1)2.11)1 corresponds to a standard Feynman diagrammatic technique with three 
bare propagators (lines in the diagrams): (vv) , given by 12.7)) . 1)2.8)1 . and the propagators of the 
scalar fields (in the frequency-momentum and time-momentum representations): 

((f(p')o = {(f'(p)o = 

= {-icr ^ + e(k)}" 1 «-» e{t-t')o-Q l cxp{-e(k)(i-t')/<7o}, (2.15) 
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where e(k) = + uok 2 + r and 9{. . .) is the Heaviside step function, and 

( W )o = 2a {^ 2 a 2 +e 2 (k)} _1 «-> -ly exp{-e(k)|t-i'|/M; ( 2 -!6) 

the propagator (ip'(p')o vanishes identically for any field theory of the type l|2.11[) . The model 
also involves two types of vertices corresponding to the interaction terms ip'(p and ip'(vdu)(p. The 
corresponding coupling constants ( "charges" ) go and Wq defined are introduced by the relations 

1/2 

A = u Q ' .go, D = wouo/ao, (2-17) 

so that by dimension go ~ l~ e and wo ~ where is £ has the order of the smallest length 
scale of our problem. More precisely, these two lengths are rather different: the coupling go in the 
Landau-Ginzburg model is conventionally related to the molecular length, while Wo corresponds to 
the Kolmogorov (dissipation) scale of turbulence. However, in the following we will be interested 
in the behaviour of the correlation functions at distances much larger than the both these lengths, 
which allows us not to distinguish them. Thus we can write 

3o ~A £ , Wo ~A«, (2.18) 

where A sets the characteristic UV momentum scale. 

By rescaling the fields, the coupling constant wo can be placed in front of the interaction term 
<p'(vd)ip in the action l|2.11|l . which is more familiar for the field theory. We do not do it, however, 
in order not to spoil the natural form of the covariant derivative, and assign the factor wo to the 
propagator (vv)o- 

3. Canonical dimensions and renormalization 

It is well known that the analysis of UV divergences is based on the analysis of canonical dimensions 
("power counting"); see e.g. [21 El- General dynamic models of the type l|2.11[l . in contrast 
to static models (like e.g. I]2.3[l ). have two scales: canonical dimension of some quantity F 
(a field or a parameter in the action functional) is completely characterized by two numbers, 
the frequency dimension dp and the momentum dimension d k F . They are determined such that 
[F] ~ [T]~ d *~[L]~ dp , where L is the length scale and T is the time scale; see e.g. Chap. 5 in book 
H3- Our strongly anisotropic model, however, has two independent momentum scales, related to 
the directions perpendicular and parallel to the vector n, and a more detailed specification of 
the canonical dimensions is necessary. Namely, one has to introduce two independent momentum 
canonical dimensions dp and d F so that 

where L± and Ln are (independent) length scales in the corresponding subspaces. The dimensions 

are found from the obvious normalization conditions djr ± = —d^ x = 1, d\ ± = -4 ± = 0, d^ ± = 
d% = 0, d% = — dt = 1, and so on, and from the requirement that each term of the action functional 
H2.11f> be dimensionless (with respect to all the three independent dimensions separately). The 
original momentum dimension can be found from the relation d k F — d F + d F . Then, based on d k F 
and d F , one can introduce the total canonical dimension dp — dp + 2d F — d F + d F + 2d F (in the 
free theory, d t oc d\ oc <9 2 ), which plays in the theory of renormalization of dynamic models the 
same role as the conventional (momentum) dimension does in static problems; cf. Chap. 5 in book 

m 

The full set of independent canonical dimensions is needed, in particular, to identify the 
completely dimensionless parameters, which only can appear as arguments in the renormalization 
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Table 1. Canonical dimensions of the fields and parameters in the model 12.111 
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constants and RG functions. Of course, existence of several independent spatial scales is not too 
exotic; it was encountered in a number of models: ferroelectrics [33] (see also section 1.17 of book 
continuous models of self-organized criticality anisotropic versions of the Kardar-Parisi- 
Zhang model [Hlj, m- axial Lifshits points [SHI and growing surfaces, driven by obliquely incident 
particle beams pTT] . 

The canonical dimensions for the model H2.11|) are summarized in table ^ including 
renormalized parameters, which will be introduced later on. From table ^ or i equivalently, from 
the relations l|2.18|l , it follows that the model is logarithmic (the coupling constants go and wo are 
simultaneously dimensionless) at d = 4 and £ = 0, so that the UV divergences in the correlation 
functions manifest themselves as poles in e = 4 — d, £ and their linear combinations or, in general, 
as singularities at e and £ — > 0. 

The total canonical dimension of an arbitrary 1-irreducible Green function T = ($ • • • <3?) i i r 

is given by the relation 

d r = d + 2 - ^JV # d*, ^ iV$d$ = N v ,d v , + N v d 9 + iV v d v . (3.1) 

Here = {N v , N v ', N v } are the numbers of corresponding fields entering into the function T, 
and the summation over all types of the fields in H3.1(l and analogous formulas below is always 
implied. 

The total dimension dr in logarithmic theory (that is, at e = £ = 0) is the formal index of the 
UV divergence Sr = dr|e=£=o- Superficial UV divergences, whose removal requires counterterms, 
can be present only in those functions T for which 5r is a non-negative integer. The counterterms 
are local, that is, in the frequency-momentum representation the counterterm to a given function 
F is a polynomial in u>, h± and k||. Since the parameters uq and ao are dimensionless with respect 
to the total dimension dp, the index Sr gives the degree of that polynomial (with the assumption 
that lu ~ k\ ~ fcy); detailed structure of the counterterms and their dependence on ao and uq is 

directly found from the corresponding partial dimensions dp and dp'". From tableland (|3.1() we 
find 

5 r = 6-3N v ,-N v -N v . (3.2) 

Dimensional considerations should be augmented by the observation that all the 1-irreducible 
functions without the field (p' (in particular, all functions involving only velocity fields) contain 
closed circuits of retarded propagators ((p'tp)o, vanish and do not require counterterms; see e.g. [^j- 
The action (f2~TT|) is even with respect to the reflection ip' — > —ip', p — > —p>, so that all correlation 
functions with odd total number of the fields tp' and tp also vanish (no diagrams for such functions 
can be constructed). It is therefore sufficient to consider only 1-irreducible functions with N v i > 1 
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and even sum AL/+i\L,. Straightforward analysis of the expression (|3.2|) then shows that superficial 
UV divergences can be present only in the following 1-irreducible functions: 

(<p'(p') (5 = 0) with the counterterm </?V'> 

(tp'tp) (S = 2) with the counterterms ip' dtp, (p'd 2 (p, (p'dj_tp, Totp'p, 
(<p'(p 3 } (S — 0) with the counterterm p'p 3 , 
(cp'tpv) (5 = 1), 

for which the counterterm necessarily reduces to the form (p'(vidi)(p = p'vd\\p. All such terms are 
present in the action l|2.11|l . so that our model appears multiplicatively renormalizable. 

The superficial divergence in the function (p'pvv) with 5 — and the counterterm (p'tpv 2 , 
allowed by the dimension, is in fact forbidden by the Galilean symmetry. Furthermore, the latter 
requires that the counterterms <p'd t <p and tp'(vidi)p enter the renormalized action only in the form 
of Lagrangian derivative tp'Vttf- 

The arguments based on the Galilean symmetry are usually applied to the velocity field 
governed by the Navier-Stokes equation, and generally become invalid for synthetic Gaussian 
velocity ensembles. It turns out, however, that for a Gaussian ensemble with vanishing correlation 
time the Galilean symmetry takes place; see e.g. |2T]. This issue, along with the consequences 
of the Galilean invariance for the renormalization in our model, are discussed in |Appcndix A| in 
detail. 

We conclude that the renormalized action can be written in the form 
S R {p',p) = S v {v) + Z x o-{<p'f + 

-Z2<jVt ( P + Z^d\}p + Ziud 2 <p — Z$T(p — Zqcju 1 / 2 [i E tp 3 /6 . (3-3) 

Here a, r, u, w and g are renormalized analogs of the bare parameters (with the subscripts "0") and 
fi is the reference mass scale (additional arbitrary parameter of the renormalized theory). Since 
the first term iS,j(v) is not renormalized, the amplitude Do is expressed in renormalized parameters 
as 

Dq = wouo/ao = wu^/a. (3-4) 

Expression is equivalent to the multiplicative renormalization of the fields ip — > <pZ v , 

ip' — * ip' Ztpi and the parameters: 

<T = &Z„, t = tZ t , it = uZ u , g a = gn E Z g , w Q = wpfiZ w (3.5) 

(no renormalization of the velocity field is needed: Z v = 1). The constants in Eqs. (|3.3J) and (|3.5|l 
are related as follows: 

Z\ = ZqZ 2 ^,, Z-2 = Z a Z {p i Z^p, Z% = Z^iZip, 

Zi = Z u Z^i Zy, Z^ = Z T Z v i Z^p, Z§ = Z g Z\j 2 Z^i Z 3 ^, (3-6) 

and from the relation Ij3.4|l one obtains: 

Z U Z W Z- X = 1. (3.7) 

The renormalization constants capture all the divergences at e,£ — > 0, so that the correlation 
functions of the renormalized model Ij3.3|l have finite limits for e, £ = when expressed in 
renormalized parameters cr, \x and so on. In practical calculations, we will use the minimal 
subtraction (MS) scheme, in which the renormalization constants have the forms Z{ = 1+ only 
singularities in e and £, with the coefficients depending on the two completely dimensionless 
parameters — renormalized coupling constants g and w. 



-ip 
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4. RG functions and RG equations 

Let us recall an elementary derivation of the RG equations; detailed discussion can be found 
in monographs CI] ■ The RG equations are written for the renormalized correlation functions 
Gr = which differ from the original (unrenormalized) ones G = ($•••$) only by 

normalization and choice of parameters, and therefore can equally be used for analyzing the critical 
behaviour. The relation <Sr($, e, /x) = £(<&, eo) between the functional (|2.11() and (|3.3() results in 
the relations 



G{e Q ,...) = Z^z"?'G R {e^,...) (4.1) 

between the correlation functions. Here, as usual, N v and are the numbers of corresponding 
fields entering into T (we recall that in our model Z v — 1); cq — {do, tq, uq, u>o, go} is the full set 
of bare parameters and e = {<r, r, u, w, g} are their renormalized counterparts; the dots stand for 
the other arguments (times, coordinates, momenta etc). 

We use T>^ to denote the differential operation [id^ for fixed eo and operate on both sides of 
the equation l|4.1|) with it. This gives the basic RG differential equation: 

{Drg + N^ + N^} G r (e, //,... )=0, (4.2) 
where I?_rg is the operation T>^ expressed in the renormalized variables: 

T^RG = £V + Pgd g + (3 w d w - j u T> u - la'Da - It^t- (4.3) 
Here we have written T> x = xd x for any variable x, and the anomalous dimensions 7 are defined as 

7i? = T>n InZp for any quantity F, (4-4) 
and the (3 functions for the two dimcnsionless couplings g and w are 

Pg = V^g = g [-e- j g ], (3 W = V^w = w[-^-^y w ], (4.5) 

where the second equalities come from the definitions and the relations 1)3.5(1 . 

Equations 1)3. 6J) result in the following relations between the anomalous dimensions 

7i = 7^ + 2 7 V ', 72 = 7<r + 7 V , 13=1 V '+1 V , 

74 = 7^+73 , 75 = 7r+73 , 76 = 7s + 7«/ 2 + 7y' + 3 7v , (4-6) 

while from (|3.7|l one obtains 

In + 7m - 7<x = 0. (4.7) 

The dimensions 71-76 are calculated from the corresponding renormalization constants using the 
definition l|4.4|l . while the RG functions entering equation l|4.3|l are easily found from the relations 
gill) and ||4~7|) : 

27^' = 7i - 72 + 73, = 73 - 71 + 72, 

7m = 74 - 73 , 7r = 75 - 73 , 7<r = 72 - 73, 

lw = 72 - 74, Ig = 71 - 72 - 373/2 - 74/2 + 7 6 . (4.8) 

The diagrams needed for our second-order calculation of the critical dimensions are presented 
in the appendix B. One can see that the leading contributions to different renormalization constants 
(and hence to the corresponding anomalous dimensions) are of different order: Z12.3 = 1 + 0(g 2 ), 
Zb,s = 1 + 0(g), Z4 = 1 + 0(w). Practical calculation of the renormalization constants and 
anomalous dimensions are discussed in |Appendix B| and here we only present the leading-order 
results for the dimensions l|4.6|l : 

7i -72 = bmg 2 + 0(g 3 ), 73 = + O(f), 76 = -3n 2 g + 0(g 2 ), 

l4 = w + Kl g 2 /6 + 0{g 3 ), 75 = - Kl g + Q(g 2 ), (4.9) 
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where we have denoted g = g/(167r 2 ), w = w/(4ir 2 ) and b = ln(4/3) « 0.287683; in counting 
the orders it is assumed that w — 0(g). For generality, we give the results for the O(N)- 
symmetric model with an iV-component order parameter in 1)2.3)) ; the additional symmetry factors 
are Ki — (N + 2)/3 and K2 — (N + 8)/9. In what follows, we will only give the results for TV = 1 
and denote the new couplings g, w simply by g, w. Then expressions 1)4.9)1 take on the form 

7i-72 = fog 2 + 0( 5 3 ), 73 = 5 2 /6 + 0( ff 3 ), l6 = -3g + 0(g 2 ), 

74- W + 5 2 /6 + 0( ff 3 ), lb = -g + 0{g 2 ). (4.10) 

5. Fixed points and scaling regimes 

It is well known that possible large-scale scaling regimes of a renormalizable model are associated 
with IR attractive fixed points of the corresponding RG equations. In our model, the coordinates 
g*, u>* of the fixed points are found from the equations 

g (g*,w*) = 0, Pw(9*,w*) = 0, (5-1) 

with the (3 functions given in (|4.5() . The type of a fixed point is determined by the matrix 

= {fiy = dfo/dgj}, (5.2) 

where (3i denotes the full set of the (3 functions and gj = {g, w} is the full set of couplings. For IR 
stable fixed points the matrix fl is positive, i.e., the real parts of all its eigenvalues are positive. 

From the definitions 1)4. 5)1 . relations (|4.8(l and explicit expressions (|4. for the anomalous 
dimensions we derive the following leading-order expressions for the j3 functions: 

P g = g [-e + 3.g + w/2] , p w = w [-£ + w] (5.3) 

with the corrections in the square brackets of order 0(g 2 ) and higher. From Eqs. 1)5.1)) and 1)5.3)) 
we can identify four different fixed points; the matrix Q appears triangular for all of them, so that 
its eigenvalues are simply given by the diagonal elements £l g = d(3 g /dg and Q w = dp w /dw: 

1. Gaussian (free) fixed point: g* = w* — 0; ft g — —e, Q w = — £. 

2. w* = (exact result to all orders), g» = e/3; il g = e, Cl w = In this regime, effects 
of the velocity field are irrelevant, the isotropy violated by the velocity ensemble is restored and 
the leading terms of the IR behaviour coincide exactly with those of the equilibrium model A. In 
particular, the basic critical dimensions do not depend on £ and coincide to all orders in e with the 
well-known static exponents rj, v for the Landau-Ginzburg model 1)2.3(1 and the dynamic exponent 
z for the model A (see e.g. pHE])- However, corrections to the leading-order asymptotic expressions 
will be anisotropic and different from those for the model A; in particular, the dependence on £ 
will appear e.g. due to the correction exponent £l w . 

3. = (exact result to all orders), w* = £; Q g — (— e + £/2), Q w — £. In this regime, the 
nonlinearity (p 3 in the stochastic equation 12.1)1 becomes irrelevant, and we arrive at the model 
of a linear convection-diffusion equation for a passive scalar field ip. For the strongly anisotropic 
Gaussian velocity ensembles of the type 12.7)) . 1)2.81) such models were investigated in detail in 
Refs. E3I13E11 (mostly for d = 2, but beyond the scope of any perturbation theory) . 

4. g* = (e — £/2)/3, w* = £; fl g = (e — £/2), fl w = £. This fixed point corresponds to a 
new nontrivial IR scaling regime, in which the both nonlinearities in the stochastic equation for ip 
are important; the corresponding critical dimensions depend essentially on the both RG expansion 
parameters e and £ and are calculated as double series in these parameters; see section [B] This 




FP 1: Gaussian 



FP 2: Model A 



Figure 1. Regions of stability of the fixed points in the model 12,111 , 



behaviour reveals strong anisotropy and belongs to a new, completely non-equilibrium, universality 
class in the sense that the equal-time correlation functions are not given by a Gibbs measure 
exp{— H(ip)} with a Hamiltonian of the type l|2.3l) . 

In figure ^ we show the regions of IR stability for all these fixed points in the e— £ plane, that 
is, the regions in which the eigenvalues VL 9tW for a given fixed point are both positive. 

In the leading-order approximation (|5.3() . all the boundaries of the regions of stability are 
given by straight lines; there are neither gaps nor overlaps between the different regions. However, 
experience with analogous two-parameter models (e.g. double expansion for the stochastic Navier- 
Stokes equation near two dimensions |39| ) suggests that such behaviour can rather be an artifact 
of the leading-order approximation: the boundaries become curved and overlaps can appear if the 
higher-order corrections in the /3 functions are taken into account. In our model this definitely 
happens for the boundary between the regions of stability of the fixed points 2 and 4, as can be 
argued without practical calculation of the corrections to the functions (|5.3[) : see section 

One can see that the interval of the most realistic values of these parameters, e = 1 (d = 3) 
and < £ < 2 (see the remark above equation l|2.9|l % ). belongs completely to the region of stability 
of the most nontrivial fixed point 4. It is also worth noting that, for all fixed points, the coordinates 
<?*, w* are positive in the regions of their IR stability, in agreement with the physical meaning of 
the these parameters: w enters the amplitude in a pair correlation function and g > is required 
for the stability of the static model l|2.3|l . 
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6. Critical scaling and critical dimensions 

Let F be some function of n independent arguments {xi, . . . ,x n } satisfying the following scaling 
relation 

F(X ai xi, . . . , X a "x n ) = X a " F(xx, ...,x n ) (6.1) 

with a set of constant coefficients (scaling dimensions) {«i, . . . , a n , cxf} and any positive parameter 
A > 0. Differentiating 16.1( 1 with respect to A and setting A = 1 gives the first-order differential 
equation 

n 

ajDj F(x x , x n ) = a F F(x x , ...,x n ), V l = x i d/dx i (6.2) 

whose general solution has the form 

F{x x ,x 2 ,...,x n ) = xl F/ai F -^J^J > ( 6 ' 3 ) 

where F is an arbitrary function of (n — 1) arguments. Obviously, the dimensions are determined 
up to an overall constant factor (replace A — > A a in 1(6.1(1 or multiply (|6.2|l by a) ; this arbitrariness 
can be fixed e.g. by setting ai = 1. If on = for some Xi, this variable is not dilated in 1(6.1(1 and 
the corresponding derivative in 1(6. 2|) is absent. 

It is well known that the leading term of the large-scale asymptotic behaviour of a 
(renormalized) correlation function satisfies the RG equation 1(4.2(1 in which the renormalized 
couplings are replaced with their fixed-point values. In our case this gives 

{D^ - 7 *£> u - 1 * a V a - 7 ;p T + N^} G m = 0, (6.4) 

where 7* = j u (g = g* , w = w* ) and so on, and Gjv# is the renormalized correlation function with 
iV$ = {N^i, N v , N v } fields. The summation over all types of fields in the last term of (j6.4|l and 
analogous expressions below is implied; cf. equation H3.1(l . 

Canonical scale invariance of the function Gjv* with respect to the three independent 
canonical dimensions (sec section |3J) can be expressed by the differential equations of the form 
(for definitcness, we consider the momentum -frequency representation) 

{V u -V a - N v i - N v } G N * = 0, (6.5) 

{V ± +V fl + 2V a + 2V T + 2V U - {N v , + N v )(d - 3)/2} G Ns = 0, (6.6) 

{By - 2V U - (N v , + N v ) 12 + N v } G Ni , - 0, (6.7) 

where T>« = knd/dk», T>± = k±d/dk±, and the data from table^have been used. 

Clearly, equation (|6.4I) corresponds to the scaling behaviour (|6.1|l of the function Gjy # upon 
the dilatation of the parameters cr, /i, u and r and at fixed momentum and frequency variables; 
equation (|6.5|) deals with dilatation of the frequencies and other parameters at fixed momenta, and 
so on. We are interested in the critical scaling behaviour, that is, behaviour of the type H6.lt in 
which all the IR relevant parameters (momenta/coordinates, frequencies/times, deviation of the 
temperature from its critical value t oc (T — T c )) are dilated, while the IR irrelevant parameters 
(those which remain finite at the fixed point: cr, /1 and u) are fixed [2 E]- Thus we combine 
the equations l|6.4|l - H6.7|l so that the derivatives with respect to the IR irrelevant parameters are 
eliminated; this gives the desired equation which describes the critical scaling behaviour: 

{V ± + A||P|| + A W P W + A r P T - JV*A*} Giv* = 0. (6.8) 

Here Aj^ = 1 is the normalization condition, while the critical dimensions of any other IR relevant 
parameter F is given by the general expression 

A F = 4+A,|4+A a ,d^+7^ (6.9) 
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with the canonical dimensions from table ^ and 

A w = 2 + 7 ;, A|| = (2 + 7 :)/2. (6.10) 

Below we will concentrate on the two nontrivial fixed points 2 and 4; see sectional 
In the fixed point 2, where u>* = exactly, all the anomalous dimensions 7^ in the RG 
equation 14. 2|) . I|4.3[) coincide with their counterparts for the equilibrium model A. In particular, 
7m — 74 ~~ 73 — an( l therefore A|| = 1; the SO(d) symmetry violated by the velocity ensemble is 
restored. Furthermore, from the fluctuation-dissipation theorem it follows that 7* = 7J; see e.g. 
[31 E] . The standard notation for this equilibrium case is 

A v = d/2 - 1 + 77/2, A T = l/v, A w = z, (6.11) 

while for <ff the aforementioned relations for jp give A v > — A v + z = d/2 — 1 + z + n/2. The 
exponents 77 and v can be found directly from the static model (12.31) : they are well known from 
the 4 — e, 2 4- e and 1/W-expansions, real-space RG (all augmented by various summations), 
high-temperature expansions for the Ising model (considered the most reliable), Monte-Carlo 
simulations. The values recommended by [21 IS] are V = 0.0375 ± 0.0025 and v — 0.6310 ± 0.0015 
(Borel summation of 5-order results). For z only two terms of the 4 — e expansion are known: 
z = 2 + 0.726(1 — 0.1885e)?7 @U|; there are also four-loop results in the real-space RG @T and 
leading-order results in 2 + £ and l/7V-expansions; see the discussion in [3] |3] and references 
therein. 

Let us turn to the fixed point 4. We are going to find the critical dimensions entering the 
equations l|6.8|l and (|6.10(l to the second order of the generalized e expansion, that is, the double 
expansion in e and £ with the convention that £ = 0(e). From the relations l|4.8[) and (|4.10(l it 
follows that J<r,tp,tp> = 0(g 2 ) and 7„ = 0(w), so that in order to find 7* / with the accuracy of 
0(e 2 ) it is sufficient to calculate the coordinate g* only to the leading order 0(e). At first sight, 
the second-order calculation of the coordinates g», u>* is needed to find the 0(e 2 ) contribution in 
the dimension 7*. However, this calculation can be avoided with the aid of the exact identity 

7„ = ? + 7:, (6.12) 

which follows from the relations f3 w = w [— £ — 7^] in Ij4.5|l and 7« + 7™ = 7<r in l|4.7J) for any fixed 
point at which f3 w = and w* 7^ 0. Then from (|6.10|) it follows that 

2A,|=A w +e (6.13) 

exactly, and from the explicit expressions (|4.10|) one obtains (e = e — £/2): 

A v = 1 -e/2 + (66+ l)e 2 /486, 

A v , = 3 - s/2 + (106 - l)e 2 /72, 

A,, = l + £/2 + (66-l)e 2 /108, 

A w = 2 + (66-l)e 2 /54, 

A T =2-e/3 (6.14) 

with corrections of order 0(e 2 ) for A T and 0(e 3 ) for the other dimensions; we recall that 

b = ln(4/3) w 0.287683. 

For the velocity field, relations 16.9|l . I|6.10JI and (|6.12(l can be combined to give the exact 
expression 

A„ = 1 + 7* - 7 ;/2 = (A w - 0/2, (6.15) 

in agreement with the explicit factorized form of the velocity correlation function (|2.7II . (|2.8|l . 

It remains to note that the critical dimensions l|6.14|l coincide up to the order 0(e) with their 
counterparts in (|6.11() at the ray e > 0, £ = 0, the boundary between the regions of stability of 
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the corresponding fixed points 2 and 4 (determined in section [3] in the first-order approximation 
for the (3 functions), but differ in order 0(e 2 ). This is a clear indication that a straight boundary 
without gaps and overlaps is an artifact of the first-order approximation: the boundaries become 
curved and overlaps or gaps appear when the higher-order corrections to the (3 functions are taken 
into account, as it happens in the analogous double expansion for the stochastic Navier-Stokes 
equation near two dimensions 39 . 

7. Discussion and conclusion 

We have studied effects of turbulent mixing and stirring on the critical behaviour of a fluid system 
(binary mixture, nematic liquid crystal) with a purely relaxational dynamics of a non-conserved 
order parameter, known as model A The velocity was modelled by Gaussian statistics with 

vanishing correlation time and strongly anisotropic correlation function cx 5(t — t')/k d ^ 1+ ^; see 
equations i|2.7J) . 1J2.8J) . Such ensembles were employed earlier in |231 1241 125] in the analysis of the 
two-dimensional passive turbulent advection (linear equation for the scalar field) . 

The model, originally described by a stochastic differential equation H2.1JI - H2.4J) . can be 
reformulated as a multiplicatively renormalizable field theory IJ2.11J) . which allows one to apply 
the field theoretic RG to study its critical behaviour. The model reveals four different IR scaling 
regimes, related with the four different fixed points of the RG equations; their regions of stability in 
the e-£ plane are identified in the leading order. These regimes correspond to: (1) Gaussian (free) 
model, (2) equilibrium critical dynamics (standard universality class of the model A, interaction 
with the velocity field is irrelevant), (3) linear passive scalar advection (the <p A term in the Landau- 
Ginzburg Hamiltonian is irrelevant) and (4) the most nontrivial strongly anisotropic scaling regime 
in which the both interactions are important; it corresponds to a new non-equilibrium universality 
class. 

It was shown that the equilibrium critical regime (model A) becomes unstable for the realistic 
range of parameters d < 3 and < £ < 2, which includes Kolmogorov spectrum (£ = 4/3) and 
Batchelor limit (£ = 2). It is replaced with the new non-equilibrium regime; the corresponding 
critical exponents are calculated to second order of the corresponding RG expansion, which in this 
case takes on the form of the double expansion in e and £; explicit expressions are given in IJ6.14J1 . 

Let us discuss the consequences of the general scaling relations, derived in section [HI for the 
most interesting special case of the pair correlation function. They result in the scaling expression 

<p(x + r, t + t') p(x, 0) = r- 2A - T (to r^ , t/r£> , HI /^" ) , (7-1) 

where r±_ = |rj_|, r\\ = \r\\ \ and J- is some scaling function. It is usually assumed that T has finite 
limits for to oc (T — T c ) = (that is, exactly at the critical point) and/or for t = (equal-time 
correlation function). Then from (J7.1|) one obtains 

<p(x + r, t) ^(x, £)) = r- 2A * f (r,| /rj 11 ) (7.2) 

with another nontrivial function T(x) = ^-"(0,0, x). The two last arguments in the scaling 
representation IJ7.1JI can also be chosen in the form r±/L±(t) and r\\/L\\{t) with two different 
characteristic length scales 

L±{t)~t a ^, L||(t)~i Q ", a± = l/A u , a||= An/A*,, (7.3) 

with the exact relation 2au = 1 +£a±_ following from equation IJ6.13JI . For the most realistic values 
e = 1 (d = 3) and £ = 4/3 (Kolmogorov spectrum of the velocity) explicit results IJtj.l4|) give 

A w w 2.0015, a ± w 0.4996 and a\\ » 0.833, (7.4) 
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while for e = 1 and £ = 2 (Batchelor limit, smooth velocity field) in the same approximation one 
obtains 

A w = 2, a_L = 0.5 and a\\ = 1, (7.5) 

with possible corrections from the 0(e 3 ) terms in (|6.14|) . It is worth noting that the 0(e 2 ) 
contributions to these results are almost negligible, so that A w appears almost indistinguishable 
from its canonical value A w = 2. On the contrary, the analog of Fisher's exponent rj = 
£/2 + (66 + l)e 2 /243, determined from tfFTty using the "equilibrium" relation (|6.11|) . markedly 
deviates from its canonical (vanishing) value due to the 0(£) term: r\ w 2/3 for £ = 4/3 and 77 w 1 
for £ = 2. This is reminiscent of the observation made in Refs. |SJE] (however, for a conserved 
order parameter and a non-random velocity) that the critical fluctuations are suppressed by the 
flow and the behaviour of the system becomes close to the mean-field limit in a strong shear; see 
also discussion in [T3] . 

Existence of two different length scales (|7.3|) with power-law dependence on the time was 
established in a number of studies within numerical simulations |15l I16j , approximate analytical 
solutions J2j and exactly soluble simplified models ^B]- As a rule, those authors dealt with binary 
mixtures in the coexistence (two-phase) region (finite and negative to oc (T — T c )) in the presence 
of a uniform laminar shear flow, while our results refer to a system near its critical point and in a 
chaotic velocity ensemble. [For finite r < 0, phase separation occurs at length scales comparable 
to the typical size of turbulent eddies, the situation which is much more difficult to achieve in 
practice in the vicinity of the critical point (to ~ 0), at least for binary mixtures; see also the 
discussion in Ref. jS].] Thus a priori one should not have expected a good quantitative agreement 
for the exponents in (|7.3|) . Surprisingly enough, our answers (I7.4f) and l|7.5|l for the exponents 
appear not inconsistent with the results a± = 0.5 and otu = 3/2, derived earlier in Refs. |17j for a 
non-conserved order parameter in a uniform non-random shear within the so-called Ohta-Jasnow- 
Kawasaki approximation 021- Although the values of an are rather different in JJj and l|7.4|l . 
H7.5(l . they are always markedly larger than a±. The same inequality a» > a± for the exponents 
was also established in two dimensions JJj and for exactly soluble models |18j . 

Let us briefly discuss the general case (|2.5|l with u^O. Nonvanishing mean velocity u gives 
rise to the additional term o~Qtp' (uidi)ip — o~oip' (ud\\)(p with u — |u| in the action (|2.11() . which 
simply results in the replacement w — > 10 — u ■ k — oj — uk\\ in the propagators H2.15JI . I|2.16(l 
and r|| — > r|| + ut in the final scaling expressions like (f T . 1 1> - This fact can be compared with the 
observation made in Refs. |16l ITTj that, for u ^ 0, the proper scaling variables are not simply 
related to parallel and perpendicular directions. The dependence on u disappears at t = 0, that 
is, in the equal-time correlation function l|7.2[l . 

Another interesting quantity is the "crossover exponent" x m the relation ST C tx Re x between 
the Reynolds number Re and the shift ST C of the critical temperature due to the mixing, with 
experimental estimates x ~ 1.4-2.1 ^2]- In the RG framework, this exponent can be identified 
U\ El as x = H^minl, where v ~ 0.63 is the classical critical exponent (|6.11f) for the Landau- 
Ginzburg model (|2.3I) and f2 m i n is the minimal (maximal by the modulus) negative eigenvalue of 
the il matrix (|5.2(l at the equilibrium scaling regime (model A or, in our notation, fixed point 2). 
In our case f2 m i n = il w = — £ for the point 2; see sectional This gives x — 1-2 for the Kolmogorov 
spectrum (£ = 4/3) and x — 1-26 for the Batchelor limit, which is better than the estimate x — 0.8 
obtained in jlOj for strongly anisotropic non-random shear but worse than the RG result x — 1-74 
obtained in for a random isotropic velocity ensemble with the velocity spectrum cx 1/fc 2 . Of 
course, the disagreement can be explained by the non-conservation of the order parameter in our 
model. [In this connection it should be mentioned that the discussion of section VI in |H] for the 
general exponent in the velocity correlation function, denoted as l/fc 2+oe in equation (6.1) of 8 , 
contains an error: the second (3 function in (6.8) must be j3\ = —\r {(1 + a)e — ...}. Thus the 
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conclusions made in the following discussion about the independence of the critical exponents on 
a (in our notation, on the relation between the two RG expansion parameters e and £) can be 
erroneous and must be revisited.] 

It remains to note that for d < 3 and not too small £, the fixed point 4 becomes unstable while 
the point 3 becomes IR attractive (see figure QJ, the ip 4 interaction in l|2.3[l becomes irrelevant, 
and the IR behaviour of the model coincides with that of the linear passive scalar advected by 
the anisotropic Gaussian velocity ensemble 12.7(1 . 1)2. S[l . For d — 2, this regime was investigated in 
detail in Refs. [2S1 EH ESI 

We may conclude that our simplified model of a non-conserved order parameter and Gaussian 
velocity ensemble captures important characteristics of a real second-order phase transition in a 
stirred fluid system: persistence of a critical scaling regime; emergence of a new non-equilibrium 
universality class with a new set of critical exponents, rather different from the classical ones; 
existence (for a strongly anisotropic velocity ensemble) of two different length scales (with a power- 
law time dependence), and so on. Further investigation should take into account conservation of 
the order parameter and its interaction with other thermodynamical degrees of freedom (mode- 
mode coupling), compressibility, non-Gaussian character and finite correlation time of the velocity 
field, and so on. This work is now in progress. 
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Appendix A. Consequences of the Galilean symmetry 

In this Appendix we will explore consequences of the Galilean symmetry for the renormalization 
of the model l(2.11|l . The models with synthetic Gaussian velocity ensembles rule, not 

invariant with respect to the Galilean transformations. Nevertheless, if the velocity is not correlated 
in time, nontrivial parts of the 1-irreducible correlation functions appear invariant (more precisely, 
see below), and the Galilean symmetry can be used to restrict the form of the counterterms. In 
this sense, the symmetry of the counterterms is higher than the symmetry of the action functional. 

For the most of the following discussion, precise form of the nonlinearity in (|2.1() is unessential; 
it is only important that it is consistent with Galilean symmetry. The velocity field will be taken 
divergence-free, Gaussian, with zero mean and the correlator 



Our model i|2.7|) , l|2.8(l corresponds to a special choice of the function , but in what follows its 
precise form is also unessential. 

Field theoretic formulation H2.11|l means that the generating functionals of total [G(A)] and 
connected [W^(A)] correlation functions of the original stochastic problem can be represented by 
the functional integral of the form 



Here and below, we denote by $ = {ip' , <p, v} the full set of fields and by A = {A^i , A v , A v } the 
full set of sources; in the expressions like 



(v i (t,x)i; i (t / > x') = 5{t - OAj(* - x ')- 



(A.I) 




(A.2) 
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summation over all types of the fields, integration over their arguments x = {t, x} and summation 
over their vector indices are always understood. All the normalization factors are included into 
the functional differential 2?$ = Dip' Dip T>v; the normalization G(0) = 1 is implied. 
The Galilean transformation is defined as 

Vi (t, X) — > Vi (t, x) = Vi(t, x + ut) - u t 

for the velocity and 

$(*,x) $(i,x) = $(i,x + ui) (A.3) 

for the other fields; here u, the parameter of the transformation, is an arbitrary constant vector. 
For the strongly anisotropic ensemble (|2.7J) . Ij2.8|) . the vector u must be parallel to v ~ n; then the 
shift of the arguments in (jA. 3|> reduces to x\\ — ► xn + u. 

The part of the action i|2.11|l which corresponds to the stochastic problem (|2.1|) at fixed v is 
clearly invariant: one has to substitute <i> — > $ and make the change of variables x + u< — > x; the 
additional terms ip'(iiidi)(p coming from the contribution with dt and from the nonlinearity cancel 
each other in the covariant combination ip'Vtf- If the velocity were governed by the stochastic 
Navier-Stokes equation with a time-decorrelated random force, the total action would also be 
invariant (see e.g. [3]), but for our synthetic ensemble (|A. 1|) variation of the action S v (v) in (|2.13|1 
is nontrivial: 

S v {v)=S v (v)+uD- 1 v + 0(u 2 ). (A.4) 
In the detailed notation, 



lD~ x v = J dt J dx J dx' u i DT j 1 (x-x')v j (t,x'), 



where the inverse linear operation for D in (|A.ljl on the transverse subspace. 

We stress that for the validity of (|A.4|> it is crucial that the correlator (|A.lJl involves the 8 
function in time. Indeed, substitution v — > v in S v produces the term 

J dt J dx J dx' u i (t,x + u*)Z>^ x (x-x / )Wj(*,x / + uf), (A.5) 

which gives S v (v) after the change of variables x + ut — > x, x' + ut — > x' due to the fact that 
the both fields in (|A.5J| have the same time argument and therefore the argument of D _1 remains 
unchanged. For a finite correlation time, expression (|A.5J) would involve the double time integral, 
the argument of D^ 1 would be shifted by u(t — t') and the original action S v in the right-hand 
side of (|A.4|I would not be formed. 

Let us make the substitution $ — * $ in the functional integral (|A.2|) . This is just a change of 
integration variables, its Jacobian equals unity, so the integral 

G(A) = J £>$ exp{5($) + A®} (A.6) 

is in fact independent of the parameter u from (|A.3J| . In particular, this means that its first 
variation with respect to u vanishes. Let us denote by S U F($>) the linear-in-u term in the Galilean 
transformed quantity F(&). Then we have 5 u v — (ud)v — u for the velocity, 5 U & — (w9)$ for 
the other fields and 5 U S(&) = 5 u S v (v) = uD^ 1 !! for the action functional. Substituting these 
expressions into the first variation of (|A.6|I gives the identity 

J £>$ {A(ud)<I> - uA v + uD- x v) exp{5($) + A$} = (A.7) 

with implied summations over all types of fields, integrations and so on. The fields can be taken 
outside the integral in (|A.7(1 as variational derivatives with respect to the corresponding sources, 
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$ — > 5/8 A, which gives for the functional W(-*4) from (|A.2|I the following differential equation in 
variational derivatives: 

A(u 9 )^- U A v+ uD-^=0. (A.8) 

It is well known that the generating functional T(<I>) of the 1-irreducible correlation functions 
(sometimes referred to as "effective action") is obtained from VK(v4) as the functional Legendre 
transform with respect to the sources A (see e.g. [2 [3]): 

m = W(A)-A*, ^ = *. ^P = -A (A.9) 

Here the sources are (implicitly) expressed as functions of the fields for a given iy(A) using the 
second relation, while the third relation explicitly determines A in terms of <f> for a given T($). 
Substituting l)A.9|) into (|A.8J| gives the following equation for r(<f>): 

^W*-JT&. =uD -*v. (A.10) 
<5$ ov 

The left-hand side obviously represents the first variation 5 l ,r($) of the functional (| A. 9|> with 
respect to the Galilean transformation i|A.3|) of its functional arguments, while the right-hand side 
is nothing but the variation 6 U S($>) of the action (|2.11|) . It is well known that the functional (| A. 9|> 
can be represented as the sum T($) = 5(<f>) + f ($) of the action <S(<i>), which contains all the 
tree ("loopless") graphs and the terms not represented by graphs, and the nontrivial part T(<I>) 
which contains all the graphs with loops (and hence those with all possible UV divergences); see 
e.g. 0Ej. From iTOTil we conclude that the total non-invariance of r(<£>) is brought about by the 
action term, while the second contribution appears invariant: (5 u r($) = 0. 

The last relation holds for arbitrary values of the model parameters, including d, £ and the 
coupling constants (|2.18l) . Therefore it remains valid in the perturbation theory and is preserved by 
the renormalization procedure. We thus may conclude that the contribution of the counterterms 
(determined by the nontrivial term f (<&)) must also be Galilean invariant, in spite of the fact 
that the total functional T($) is not. This justifies the statements made in the analysis of the 
renormalization of our model in section |3J the counterterm (p'ipv 2 , allowed by the dimension, is 
not invariant and therefore it is forbidden; the counterterms ip'dttp and ip'(vidi)ip can appear only 
in the form of the Galilean covariant combination ip'\7 t ip. 



Appendix B. Calculation of the Feynman diagrams 

In this appendix we will briefly discuss the main points concerning the calculation of the 
renormalization constants and the corresponding Feynman diagrams. In order to find the 
anomalous dimensions (|4.6|) in the approximation (|4.9|) , (|4.10f) , one has to calculate the 1-irreducible 
correlation functions (if' if'), ((/?V) to the two-loop order and (ip'ipipip) to the one-loop order of the 
renormalized perturbation theory. The corresponding diagrammatic expressions are given in figures 
IBlllB3l fwe do not show some of the diagrams, which are topologically possible but vanish because 
of special reasons; see below). 

The wavy lines denote the pair correlator of the velocity (|2.8|) . The bare propagators l|2.15|l . 
H2.16|l are denoted by solid lines: the line without a slash denotes the propagator (iptp)o while 
the line with one slash denotes (ip(p')o; the slashed ends correspond to the response field ip' . The 
propagator (ip'Lp')o vanishes in any dynamic model of the type (|2.11|) . so the solid line with two 
slashed ends does not occur in the diagrams. There are two types of vertices (shown by thick dots): 
the quartic one corresponds to the interaction tp'ip 3 and the triple one corresponds to ip'(vd\\ )ip. 
The symmetry coefficients are shown for the one-component scalar field (N — I). 
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(ip'ipipip) = Z 6 



+ 3 




+ 



No 12 



+ 
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+ ... 



No 13 



No 14 



Figure B3. 1-irrcduciblc function (ip'iptpip): One-loop approximation. 



All the diagrammatic elements should be expressed in renormalized variables using the 
relations (|3.3[l l|3.7|l . However, in our approximation the Z's should be retained (with the 
appropriate accuracy in g and w) only in the bare terms of the functions ($...$), while in 
the diagrams they should be replaced with unities. In other words, the passage to renormalized 
variables in the diagrams is provided by the simple substitutions uq — > a, uq — > u, tq — > r, go — > g\i z 
and wo — ► wfjfi . 

In the practical calculation we use the MS scheme, where the renormalization constants are 
independent of the specific choice of the IR regularization. It is then possible to calculate the 
constants directly in the critical ("massless") model, that is, at r = in the renormalized analogs 
of the lines (|2.15() . I|2.16[) . Then, in the calculation of the constants Zi with i ^ 5, the diagrams 
involving self-contracted propagators (yy)o ca- n be treated as zero and, for this reason, are not 
shown in the figures. The only exception is made for the one-loop diagram No 7 with a self- 
contracted solid line in the function (<pip') needed for the calculation of Z5, in which r should be 
retained. Depending of the type of a diagram and specific way of calculation, the IR regularization 
is either provided by the sharp cutoff or is not needed at all (see below). The diagrams with self- 
contracted lines (tp(p')o should also be discarded according to the general rules of the diagrammatic 
technique for the dynamic models of the type 12.11JI : see 

Owing to specific properties of our model, many diagrams still shown in the figures vanish 
identically or appear UV finite and therefore do not contribute to the renormalization constants. 

Since the velocity field is transverse (divergence-free), the derivative du at the triple vertex can, 
if desired, be moved onto the response field using the integration by parts; see Eq. I|2.12|l . Thus 
in any diagram involving TV external vertices of the type <p'(vd\\)(p, the factor p^ with external 
momenta pu will be taken outside the corresponding integrals over the internal momenta and 
frequencies. This reduces the dimension of the integrand and can make the diagram UV finite. In 
particular, this makes UV finite the diagrams Nos 1, 3, 4, 13 and 14 (the latter in fact vanishes, 
see below). 

The diagrams Nos 5, 10, 11 and 14 effectively involve closed circuits of retarded propagators 
(<p<p')o (self-contracted chains of step functions; see Eq. I|2.15|l ) an d therefore also vanish. Such 
effect is well known for dynamic models of the type p.lljl and is related to causality; see [3 EJ- 
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However, for the diagrams No 11 and 14 it is crucial here that the propagator of the velocity 
involves the 5 function in time. These arguments, however, do not apply to the one-loop diagram 
No 6, which requires more accurate consideration (see below). 

There are more diagrams in the functions (ip'tp 1 ) and (ip'ip) having the same topology as the 
diagrams Nos 4 and 5, with another placements of the slashes. They also vanish because of the 
two reasons explained above, and we do not show them in the figures. 

As a result, in the functions (ip'tp 1 ) and (ip'ipipip) only the diagrams Nos 2 and 12 are divergent. 

Integration over times (or frequencies) in the diagram No 9 leads to the expression (up to a 
numerical factor and with implied IR cutoffs) 

pn [ dk [ dq S JJV3 (B.l) 

"./ J [(p + k) 2 + q 2 + (q + k) 2 ] q2k d ± - 1+i 

where one external momentum p\\ has already appeared as an overall factor. Due to the presence 
of the 6 function in the integrand, one can replace (p + k) 2 — > p?, + (pj_ +k^) 2 in the denominator. 

Thus the expansion of the integral in small momenta begins only with quadratic terms: p\ (due to 
the persisting SO(d — 1) invariance) and p|. The total expression IjB.ljl can only begin with cubic 
terms and is therefore UV finite, because the possible superficial divergence of the function (ipip 1 ) 
must be quadratic; see section|21 It remains to note that in the isotropic case the UV finiteness of 
such diagram would be guaranteed by the transversality of the velocity propagator; see the remark 
in Ref. 0. 

Thus we are left with the five UV divergent diagrams Nos 2, 6, 7, 8 and 12. 
The analytic expression for diagram No 6 has the form 

ftM = torf J* } J* j ,D.Vi ) + ' e(p _ k) (B.2) 

with D v from (|2.8[l . the prefactor coming from the vertices and e(k) = k^+uk 2 coming from l|2.15|l : 
the result is independent of the external frequency. Integration over to involves the indeterminacy 

' =a- 1 6(0), (B.3) 



(2tt) — \auj + e(p — k) 

where 0(0) is the step function at the origin. This indeterminacy reflects the details of the velocity 
statistics lost in the white-noise limit H2.7(l and it should be carefully resolved; see e.g. the discussion 
in the appendix of Ref. In our case, the function 5(t — t') should be understood as the limit 
of a narrow function which is necessarily symmetric in t t', because (12 .711 is a pair correlation 
function. Thus the quantity in 1B.3|) must be unambiguously defined by half the sum of the limits: 
#(0) = 1/2. Then after the trivial integration over k|| and using (|3.4|> one obtains 

2 wufJ,t f dk± 1 
D ^ = ~ P l\ — J (2^) • 

Finally, the integration over kj^ gives 

S 1 1 

A»(p) = -plwu^/mf 2(2 ^L) j = -pf ^ j + UV finite part, (B.5) 

where Sd = 2-K d / 2 /T(d/2) with Euler's V function is the surface area of the unit sphere in d- 
dimensional space. 

One important remark is in order here. In models with a single UV regulator (say, model A 
with e) the UV singularities manifest themselves as poles in e, and the MS scheme is defined such 
that all the renormalization constants have the form Z = 1+ only poles in e. In models with two 
regulators, like e and £ in our case, there are subtleties in defining the MS scheme: for example, 
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is the ratio e/£ a pole or a finite quantity. The final physical results must be independent of the 
choice of the renormalization scheme. Practical calculations in analogous two-parameter models 
(e.g. two-loop calculations for the stochastic Navier-Stokes equation near two dimensions |39p 
confirm that this is indeed true. In our calculations we always assumed that e ~ £, treated the 
combinations like e/£ as UV finite and did not include them into the renormalization constants; see 
the last equality l|B.5(l . However, in the leading-order approximations these subtleties are not too 
important. In particular, another (and eventually equivalent) possibility is to include the factor 
5'd_i/2(27r)^~ 1 ^ (and not only its value at e = 0) into the definition of the new coupling constant 
w\ see the text below Eq. (|4.9|l . This will include all powers of the ratio e/£ from l|B.5() into the 
corresponding renormalization constant without changing the anomalous dimension 74. 

The remaining diagrams Nos 2, 7, 8 and 12 do not involve the velocity propagator (w)o and 
can be reduced to the well-known diagrams of the isotropic model A. Consider the diagram No 8 
as an example. The corresponding analytic expression can be represented in the form 

\ 2 I{u, p\ + up\, u) = A 2 u _1 I(u;, p\+up\, u = 1). (B.6) 

Here, u> and p = p± + p|| are the external frequency and momentum, A 2 comes from the vertex 
factors and /(. . .), after the integrations over the internal times or frequencies, is represented as 
a double integral over the two integration momenta, say, k and q. All the momenta enter the 
integrand via the functions of the type e(k) = +wk| coming from the expressions (|2.15l) . (|2.1(il) 
with the replacements uq — > u and to — > (see the discussion in the beginning of the appendix). 
This explains the fact that /(. . .) depends on p through the only scalar argument p\ + up?,. This 
also gives the second equality in (|B.6|I after the rescaling of the integration momenta ukii — * k|| , 
uc l\\ ~^ 9.11 • The factor u _1 arises from the Jacobians and combines with A 2 to give the prefactor 
A 2 u _1 = g 2 ; see Eq. (|2.17|) . However, the quantity on the right-hand side of (IB.6|) is nothing other 
than the analytic expression for the diagram No 8 in the isotropic (u = 1) model A, while the 
dependence on u persists in its modified momentum argument p 2 — * p\ + up 2 . 
The isotropic integral is represented in the form 

p, u = 1) = {-io-uA +p 2 B} - + UV finite part 

with known (see e.g. 00]) dimensionless coefficients A, B, which for (|B.6|I gives 

\ 2 I(lu, p\ + upp u) = g 2 j-ioWl + (p\ + up 2 )B^ - + UV finite part. (B.7) 

This expression fully determines contribution of the diagram No 8 to the renormalization constants 
■^2,3,4 hi our model l|3.3|l . Similar considerations determine the contributions of the diagrams 2, 7 
and 12 to the constants •Z'1,5,6 m terms of the known coefficients for the model A. The latter are 
well known, but for completeness we will discuss the corresponding calculational techniques which 
proved to be useful in the three-loop calculation in the model A 0U] and might be interesting in 
itself. 

Two key points are as follows: the convolution of two functions of the form 

F (a; a) = {-iua + k 2 )- a (B.8) 

is a function of the same form, 

F{a; a) * F(f3; b) = 6{ab) K 2 (a, (3; a, b) F{a + I3- d/2 ~l;a + b) (B.9) 

with the coefficient 
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while the product of two such functions can be represented as a single integral of a function of the 
same form with the aid of the generalized Feynman formula: 

F(a; a) ■ F{/3; b) = / ds s a ~\l - sf~ l F{a + (3;as + b{\ - «)). (B.10) 



r(a)r(/3) Jo 

Equation (|B .91) can be obtained from the fact that in the {t, x} representation the function l|B.8|) 
takes on the form 

[a > ' (47r) d / 2 r(a) P \ 4tJ' [ti - U) 

and the product of such functions (which corresponds to the convolution of their Fourier transforms) 
is obviously a function of the same form. Note that for ab < the convolution l|B.9jl vanishes 



because it corresponds to the product of a retarded and an advanced functions of the form (|B.11(1 
In the same manner, for the convolution of three functions IB. 81 one obtains 



F(a; a) * F(j3; b) * F( T , c) = K 3 (a, f3, 7; a, b, c)F(a + /3 + 7 - d-2) (B.12) 
with the coefficient 

*c a 7; «, », = *«.wo r ^±i^^ x 

X0 J/2-., (> J/2-« ( .J/2- 7(a + j + I . ) » +/ J+-,-3J/2-2^ 

and so on. 

In the notation (|B.8(I . the bare propagators of the scalar fields ip, ip' can be written as 
(cp<p')o=F(l,o-), (<p'<p) = F(l,-o-), ( W ) = / duF(2;u). 



In the last expression we used the equation (|B.10|I with a = —b = a, a = = 1 and introduced 
the new integration variable as u = er(l — 2s). 

Consider again the diagram No 8 as an example of the calculation of the pole part in e. It 
is convenient to set a = 1 (the dependence on a can easily be restored by dimensionality) and 
discard the factors like (47r) _d / 2 which always combine with g to form the new coupling constant 
g; see the text below Eq. (|4.9I) . Then the analytic expression for the diagram No 8 reads (we omit 
the factor g 2 ): 

D& = J dm J du 2 F(2,u 1 )*F(2,u 2 )*F(l,l). (B.13) 
Applying the reference formula (|B . 1 2|l to the integrand of (|B . 13|> gives 
Dg = r(3 — d) / du\ I du2 u\ d ^ 2 u 2 , d ^ 2 x 

x (tti +u 2 + lf- 3d / 2 F(3 -d, Ul +u 2 + 1). (B.14) 

Note that only positive values of 1*1,2 survive in i|B.14(l due to the 9 functions in (|B.12(I . 

In order to extract the pole part of (|B.14(1 . it is sufficient to replace r(3— d) = T( — 1+e) —> —1/e 
and to set e = in the integrand; this gives: 

D s = -— [-iul 2 +p 2 h] + UV finite part, (B.15) 
with the coefficients 

r2 



In = dm du 2 ■ — — = / dup(u) — — — (B.16) 

(u 1 +u 2 + l) n J {l + u) n 
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and the weight function 

for u < 1. 



p(u) = 



for u > 1. 



The integrals in (|B.16|I are easily calculated directly, but it is instructive to use the following trick, 
which proved extremely useful in the higher-order calculations. Let us introduce the generating 
function 

f 2 1 

I(z) = / dup(u) . (B.17) 



It is easily seen that the integrals I2 and ^3 from l|B.16|) are the second and the third coefficients 
in the expansion of the quantity 1(1 ~ 5) in 5. Calculation of the integral (|B.17(1 up to irrelevant 
constant terms gives 

z(z + 2) z + 2 
l(z) ~ 2 In — - + 2 In -, 

W (z + 1) 2 2 + 1 

and the expansion 

1(1 -S)= const + (51n(4/3) + 5 2 /6 + 0(6 3 ) 



gives the desired coefficients I2 3 in (|B.15() . which in their turn determine the coefficients A, B in 

Once a renormalization constant Zi has been calculated, the corresponding anomalous 
dimension is readily found from the relation 

7i = (f3gd g + (3 w d w ) In Z t = -(eV g + £D W ) In Z t . (B.18) 

In the first equality, we used the definition (|4.4p . expression (|4.3|) for the operation T)^ in 
renormalized variables, and the fact that the Z's depend only on the two completely dimensionless 
coupling constants g and w. In the second equality, we retained only the leading-order terms in the 
P functions H4.5JI . which is sufficient for our approximation. The factors e and £ in (|B.18J| cancel 
the corresponding poles contained in In Zi, which leads to the final UV finite expressions for the 
anomalous dimensions, given in (|4.9|l and (|4.1U|) . 
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